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Abstract. In this paper, for continuous, linearly-controllable quadratic control systems with 
a single input, we propose an explicit, constructive method for studying their Brunovsky forms, 
initially studied in In this approach, the computation of Brunovsky forms and transformation 
matrices and the proof of their existence and uniqueness are carried out simultaneously. In addition, 
we show that quadratic transformations in |2| can be simplified to prevent multiplicity in Brunovsky 
forms. We further extend our method for studying discrete quadratic systems. Finally, computation 
algorithms for both continuous and discrete systems are summarized, and examples demonstrated. 
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1. Introduction. Linear control systems can be continuous in time t, 

(1.1) i^A^ + bfi, 

or discrete 



(1.2) 



where coefficients A e SR"^" and b e 3?"^^ are generally constant, and state variable, 
^(g 3fJ"), and control variable, /i(£ 3?), are continuous or discrete respectively. When 
controllable, both 1)1.1(1 and ((1.2|l admit the Brunovsky form which is derived from 
controller form, under the following linear change of coordinates and state feedback 



(1.3) 



T 
X U, 



where x and u are new state and control variables respectively, and T e 3?"^" and 
V £ 3?"^-"^ are the transformation matrices (vectors) (refer to Chapter 3 of [2]). In the 
linear Brunovsky form of Hl.l|) and H1.2|l . A and b have the following forms: 
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In an attempt to extending the Brunovsky forms for non-linear systems, which 
generally do not even have controller forms, Kang and Krener |2| studied continuous 
linearly controllable quadratic control systems with a single input, which can be writ- 
ten as (We use different notations from [3] for the purpose of clearly presenting our 
method of computation.) 



(1.5) 
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where S^^^l (^) = • ■ • , C^-fnO"^ ^ vector of n quadratic terms with symmetric 

nxn matrices F^'s (i = 1, • ■ • , n), G^fi = (Gi^/i, • • ■ , G„^^)"^ is a vector of n bihnear 
terms with G £ 3?"'^", and 0{£,,fj,)^ includes aU terms with a + b >3. Moreover, 
the foUowing additional assumptions are made for this system. First, the coefficients 
in H1.5|l and (|1.9() A, b, Fi (i = 1, • ■ • , n), G and h are assumed to be time- invariant. 
Second, the two systems are assumed to be linearly controllable; by linearly control- 
lable we mean that the linear parts of the two quadratic systems are controllable and, 
as a result, the linear parts can be transformed into the Brunovsky form with 

Kang and Krener |3] first defined quadratically state feedback equivalence up to 
second-order, or quadratic equivalence for short, under the following quadratic change 
of coordinates and state feedback, 

(1.6) C = 2: + *p[^l(x) +0(a;)^, ly = i^i + x^Qx + rx^i + 0{x, , 

in which — {x^ Pix, • • ■ , x^ Pnx)^ is a vector of n quadratic terms, and trans- 

formation matrices include symmetric Pi £ 5R"^" {i — 1, • • • , n), symmetric Q G 3?"^", 
and r G SR'^^". We can see that (|1.6|) are equivalent to 

(1.7) C = a; + «p[^l(x)-|-0(a;)^ ^ v - x^ Qx - rxv + 0{x,vf , 

and hereafter refer to transformations of 3 as (|1.7|l . Then, from all quadratically 
equivalent systems of a general system (|1.5(l . two types of Brunovsky forms were 
defined in In type I forms, the nonlinear terms are reduced into a number of 
quadratic terms ; i.e., there are no the cross terms in Xi and Xj (i ^ j) or bilinear 
terms, xv. In type II forms, only bilinear terms are kept. In both types of Brunovsky 
forms, the number of non-zero nonlinear terms is n{n— l)/2, compared to n^(n + 3)/2 
for a general quadratic system. 

In this paper, we propose a new method for studying the Brunovsky forms, first 
for (|1.5|l under transformations l|1.7|l . This method can be carried out in three steps: 
first, we find the relationships between the coefficients of H1.5|l . the coefficients of 
its quadratically equivalent systems, and the corresponding transformation matrices; 
second, from these relationships, we derive a mapping from the coefficients of (jl.5|l 
and its equivalent systems to the transformation matrix Pi, which can be considered 
as the necessary condition that all equivalent systems of Ijl.SI) should satisfy; third, we 
show how to compute, from the necessary condition, Brunovsky forms as well as the 
corresponding transformations. With this method, we find that l|1.5|l admits the same 
two types of Brunovsky forms defined in In contrast, our approach, constructive in 
nature, is capable of computing the Brunovsky forms and the transformation matrices 
simultaneously. 

We further show that quadratic transformations, (|1.7|l . can be simplified by setting 
r = to 

(1.8) C = a; + «p[^l(a;) +0(x)^ ^l = u ~ x^ Qx + 0{x,uf . 

Still defining the quadratic equivalence in the sense of 3,, the new transformations 
prevent multiple solutions of type I or type II Brunovsky forms. I.e., the Brunovsky 
form of each type and the corresponding transformation matrices Pi (i = 1, • ■ • , n) 
and Q are uniquely determined by the original system. Moreover, we apply the 
same method, but with the simplified transformations defined in (|1.8() . and study the 
following discrete system 

(1.9) e(i + 1) = Am + b^i{t) + y[2i(e(i)) + Gm^^{t) + vc*) + o(e, f^f, 
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where, similarly, ^jl^l ^-nd G£,{t)^{t) are a vector of quadratic terms and a vector 

of bilinear terms respectively. ^ We find that for p.9|) there exists only one type of 
Brunovsky form consisting of n(n+ l)/2 bilinear terms, which corresponds to type II 
Brunovsky forms of continuous systems. 

The rest of this paper is organized as follows. After reviewing the Brunovsky form 
of linear systems (Section 2), we study the Brunovsky forms of continuous quadratic 
systems in Section 3 and of discrete quadratic systems in Section 4. In Section 5, we 
summarize our method into two computation algorithms: one for continuous systems, 
and the other for discrete systems. We conclude our paper in Section 6. 

2. Review: computation of the Brunovsky form of a controllable lin- 
ear control system. We first review the computation of the Brunovsky form of a 
continuous controllable linear control system Hl.l() . ^ The procedure and the results 
also apply to the discrete system H1.2|l . Computation of the Brunovsky form for Hl.l|l 
comprises of two steps: first, the linear system is transformed into the controller form 
with a linear change of coordinates given by the first equation in (|1.3|) : second, the 
controller form is further reduced into the Brunovsky form with a state feedback given 
by the second equation in 1)1. 3|l . 

The control matrix for is defined as 



(2.1) 



C = [A''~^b 



Ab b]. 



Since 1)1.11) is controllable, rank(C) is n and, therefore, C is invertible. Denote d as 
the first row of C^^, then we can compute the transformation matrix T as 



(2.2) 



d 

dA 

dA""-^ 
dA"'^ 



and with the linear change of coordinates ^ = Tx, (|l.l)l can be transformed into the 
following controller form. 



(2.3) 
in which 



i — Ax + 



(2.4) A^T^'^AT - 







1 



1 
: 

-Vi -V2 ■ ■ ■ 



1 

-Vr, 



b = T-'b 



^Note that the discrete system 11.91 contains a term quadratic in the control variable, hp?(t), 
where h G 3?"*<l, and OJ has n^(n + 3)/2 + n non-zero nonlinear terms. 
^For further explanation, refer to 1 and Chapter 3 of 1^. 
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Then, using a linear state feedback /i 



x^v, in which, 



(2.5) 
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we can transform the controller form ^2,.'6\ into the following Brunovsky form. 



(2.6) 



From the derivation above, we can see that the Brunovsky form exists and is 
unique, and the transformations can be explicitly computed. In addition, this proce- 
dure can be easily integrated in applications related to the Brunovsky form. In the 
same spirit, we propose a method for directly computing quadratic Brunovsky forms 
and the corresponding transformations in the following sections. 

3. Computation of continuous quadratic Brunovsky forms. Note that, 
with linear transformations in quadratic systems H1.5(l and (|1.9|l will not change 

their forms except their coefficients. Without loss of generality, therefore, we hereafter 
assume that A and h in l|1.5() ( and also in (|1.9(l 'l have already been transformed into 
the forms defined by H1.4|l . In this section, we study Brunovsky forms of continuous 
system (|1.5|l under transformations (|1.7() . in a constructive manner. First, we find 
relationships between coefficients of H1.5II . coefficients of its quadratically equivalent 
systems, and the corresponding transformation matrices. Second, from these relation- 
ships, we derive a mapping from the coefficients of H1.5|l and those of the equivalent 
systems to transformation matrix Pi; this mapping is a necessary condition that all 
the equivalent systems should satisfy. Third, we show how to obtain two types of 
Brunovsky forms as well as the corresponding transformations. 

3.1 

L : - 



Quadratically equivalent system of (|1.5|) . Definition 3.1. Assume 
the linear Brunovsky form given by il.4\ ), we define a linear operator 



(3.1) 



LP 



P. 

A^P + PA, 
L ■ VP, i 



0,1,- 



Properties of L. The linear operator L has the following properties: 

1. VP = when « > 2n - 1. 

2. The nullity of L is n; if P G ker{L); i.e., LP — 0, then P can be written 

as 



P. 



0, 



in which pi , 



(-l)>,+j_„, 
,Pn are independent. 



i + j < n, 
otherwise, 
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3. L is not invertible. 

4. If P is symmetric, LP is symmetric. 

Theorem 3.2. The continuous quadratic system ^.5\) is equivalent, in the sense 
of 0', under the quadratic transformations given by j j. ?| ), to a quadratic system, 
whose ith (i = \, ■ ■ ■ ,n) equation is 



(3.2) 



Xij^i + biV + FiX + Gixv + 0{x, v)^ 



where Xn+i(t) = is a dummy state variable. Pi, ■ ■ ■ ,Fn are symmetric, Gi is the ith 
row of the matrix G, and the coefficients P^ (i — 1, ■ ■ ■ ,n) and G are defined by 



(3.3) 
(3.4) 



p = p + p+i ~LP-6.g, 



where P„+i G 3?"^" is a zero dummy transformation matrix, and bi is the ith element 
of b defined in j 

Proof. Plugging the transformations defined in (|1.7|) into H1.5f) . we obtain 

X + (x) ^ Ax + biy + ^[21 (x) + Gxz^ 

+A<p[2l(x) - te^Qx - 6ra;j/ + 0{x, vf, 

of which the ith (i = 1, • • • , 7i) differential equation can be written as 

Xi + — (x'^Pix) — Xi+i + biV + x'^FiX + Gixv 
dt 

+x'^ Pij^ix — bix'^Qx ~ birxv + 0{x, v)^ . 

After plugging x = Ax + bv + 0{x,iy)^ into the term -^{x^Pix) and collecting all 
terms whose orders are higher than two, we then obtain the equivalent system (|3.2f) 
with the coefficients given in (|3.3|l and 1)3. 4|l . ■ 
Remarks Equations 13.31 and l3.4l Dresent the relationships between the coefficients of 
H1.5|l . those of its quadratically equivalent system (|3.2(l . and the correspond- 
ing transformation matrices. We further simplify these relationships in the 
following subsections. 

3.2. A necessary condition for quadratically equivalent systems. Lemma 
3.3. The mapping from the coefficients, G and G, and r to the transformation ma- 
trices Pi (i = 1, ■ ■ ■ ,n) is given by 



(3.5) 



[P2 



[Pn 



(n) 



'-G^'-G-hr, 
2 2 2 ' 



where the operator [■](«) takes the nth row from its object. 

Proof. Since 6; = (i = 1, • • ■ , n — 1) and &„ = I, we obtain, from (|3.4|) . 



in) 



yP 



:G,; 



:G, 



-6,r. 
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Thus we have ■ 
Lemma 3.4. The mapping from the coefficients Fi and Fi (i — 1, • ■ • , to the 
transformation matrices Pi (i = 1, ■ ■ ■ ^n) and Q is given by 

(3.6) P,+i=£^Pi-^FF,_j+^FF,_j, i = l,---,n-l, 

n-l 

(3.7) Q = Y. ^'(Pn-j - K-j) - L^Pi- 

3=0 

Proof. When i — 1, ■ ■ ■ ,n — 1, from (|3.3|l , we have 

P,+i=LP,-F, + F,. 

After iterating this equation with respect to i, we obtain 13.6|l . When i = n, (|3.3|l 
can be written as 



Q Fjl Fji \jPyi. 

Since P„ is given in (|3.6|) . we then have (|3.7|l . ■ 
Definition 3.5. With A e 3?"^" and i defined in Jj.^| ) and Definition 15*. 71 
respectively, we define a series of linear operators Xi : 3?"'*" — > 5R"^" (z = 0, 1, • • -j 



Properties of X;. The hnear operators (i = 0, 1, • • •) have the following proper- 
ties: 

1. Xg transforms a diagonal matrix into a skew-triangular matrix of the 
following structure. For a fcth (fe — — n-|-l,---,n — 1) diagonal matrix 
P, we denote the fcth diagonal elements by pi = P{\k\-k)/2+i,{\k\+k)/2+i 
{I — — |fc|). Then XqP becomes a diagonal matrix with the 
following properties: (XoP)„_(|fe|_fe)/2-;+i^(|fc|+fc)/2+i = Pi, all (XoP),j 
for i+j = n + k + l + 2m, where m > 1 and n + fc + 1 + 2m < 2n^, are 
determined by pi, and the other entries are zeros. 

2. From the preceding property, Xq transforms a lower-triangular matrix 
into a full matrix and transforms an upper-triangular matrix into a lower 
skew-triangular matrix A, defined by 

(3.9) Aij = 0, when i+j<n + l. 

3. From Properties 1 and 2, we have that the nullity of Xq is 0. Therefore, 
Xq is invertible. 

4. From the definition of L, X^ = when i > n. 



(3.8) 



XoP 



[L"P] 



•^Note that, for ra = 2, there is no solution to m. In this case, all entries are zeros except (XqP) 
for i+j = n + k + l. 
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5. From the definition of Xj (i — 1, • • • , n — 1), X^P can be obtained by 
shifting XqP down by i rows and replacing the first i rows by zeros. From 
Property 1, therefore, X^ transforms a diagonal matrix P, whose main 
diagonal elements are pi , • ■ • , p„ into a lower skew-triangular matrix of 
the following structure: pi.iP)n-i+i+i,i — pi for I — i + I, ■ ■ ■ ,n, and all 
other elements except {'KiP)ij for i + j ~ n + i + 2m, where integer 
m > 1 and n + k + 1 + 2m < 2n, are zeros. 
Theorem 3.6. The mapping from Fi, G, Fi, G (i — 1, • • • , ti — 1 j, and r to Pi 
is given by 



(3.10) 



Pi 



^n-l 



1 



^ X,F, + -G 

V i=l 



- 1 - 1 ^ 

Y^X.,F,--G--br 



2 2 

Proof. From (|3.6|l we can find the nth row of P^+i (i = 1, • • ■ , n — 1) as 



[P,+i](„) = [L'Pi](„)-^[L^P,_,](„ 



i-l 



3=0 j=0 



Hence, we have 



(3.11) 



[Pl](n) 



[P2 



= xoPi -;^x,p, + ^x,p,. 



From H3.5|l and H3.11|) . we obtain 



n~l 71—1 

XoPi = ^ X,P, + :,G-Y, - -G - -br. 



Multiplying both sides by the inverse of Xq, we then have H3.10|l . ■ 
Remarks. Note that Pi is assumed to be symmetric. Therefore quadratically equiv- 
alent coefficients Fi (i = 1, • • • , rt — 1) and G have to ensure the symmetry of 
the right hand side of H3.10|l . Thus (|3.1Q|I constitutes a necessary condition 
for all equivalent systems of (|1.5|l . including the Brunovsky forms. 

3.3. Computation of the Brunovsky forms and the transformation ma- 
trices. In this section, given r as well as the coefficients of (|1.5|) . Fi {i — 1, • • • ,n) 
and G, we show how to choose Fi {i = 1, ■ ■ ■ ,n) and G in Brunovsky forms, which, 
first, satisfy the necessary condition H3.10|l and, second, have the smallest number of 
non-zero terms. 

Since P„ does not appear in the necessary condition, (|3.1U|) . we simply let P„ = 
in Brunovsky forms. To determine other coefficients, we decompose the terms related 
to original system (|1.5|) and r on the right hand side of H3.10|l as 



(3.12) 




where L, -D, U are strictly lower triangular, diagonal, and strictly upper triangular 
matrices respectively. In order to add Xq ^(^"Z^^ X^Pj -I- ^G) to the right hand side 
of H3.10|l to obtain a symmetric matrix, we have, according to the property of the 
decomposition, the following three cases. 
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1. If L — = 0, L + D + U is already symmetric, and we simply set Fi 
[i = 1, • • ■ , n — 1) and G to be 0. In this case, the Brunovsky form of (|1.5|) is 
a linear system; i.e., H1.5|) can be linearized. 

2. When L - f7^ 0, we can set Xq ^(^^^i ^iFi + ^G) to a lower-triangular 
matrix, L — U'^ . In this case. Pi — + D + U. According to the properties 
of Xq, J2^=i ^iPi + 5^ is a full matrix. 

3. When L — U"^ ^ 0, we can also set ^(^"J"^^ X^^i + if?) to be an upper- 
triangular matrix, U — L^, 'Yl^Zi '^iPi + is a triangular matrix defined 
by H3.9|l . which has n{n — l)/2 non-zero terms as U — . In this case, we 
have 

n—l ^ 

(3.13) ^XiF,-f -G = Xo(C/-L^) = Ai, 

and the first transformation matrix is given by 

(3.14) Pi=L + D + L^. 

4. In addition to the aforementioned cases, Xq ^(^"^-^^ X^fi -I- ^G) can be also 
either L — U'^ or [/ — L"^ plus an arbitrary symmetric matrix. 

Comparing Cases 2, 3, and 4, we can see that the solutions in Cases 2 and 4 have 
more non-linear terms than in Case 3. Therefore, in the Brunovsky forms, we set 
Xq ^(^"^-^^ X-iFi + iG) to be J7 — , and from H3.13|l . we can obtain two types of 
Brunovsky forms as follows. 

First, by setting G = 0, we have from (|3.13|l that J2^=i '^i^i — which is 
a lower skew-triangular matrix. To obtain the smallest number of, i.e., n(n — l)/2, 
non-zero items in the Brunovskey forms, we can select Fi to be a diagonal matrix with 
main diagonal elements as (0, • • • , 0, • • • , /»,«), which can be uniquely computed 

as follows. From properties of X^ (i = 1, • • • , n — 1), we first have 

(3-15) fi,j = (XiFi)«-i-F2j = (Ai)„_j+2j,j = 2,---,n, 

and 

(3.16) A2 = Ai-Xii?i. 

Then, we can compute f2j ior j = 3, ■ • ■ ,n from A2 and A3 = A2 — X2F2 in the 
same fashion. By repeating this process, we can obtain all non-zero elements in Fi for 
i — 1, ■ ■ ■ ,n — 1. Since all these matrices are uniquely determined by Ai, the original 
system is uniquely equivalent to the following system, 

n 

(3.17) Xi = Xi+i + biiy + ^ fijX^j +0{x,h')^, i = l,---,n, 

which is type I, complete-quadratic Brunovsky form in 

Second, by setting Fi =0 {i — 1, ■ • • , n — 1), we have G = 2Ai and the corre- 
sponding equivalent system, 

n 

(3.18) Xi = Xi+i + biU + ^ GijXjV + 0{x,iy)^, i=l,---,n, 

j=Ti-(i+l) 
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which is type II, bihnear Brunovsky form in We can see that G is also uniquely 
determined by the original system. 

Once we have the Brunovsky forms of a continuous quadratic system, we can 
solve the corresponding transformation matrices Pi (i = 1, • • • , rt) and Q from H3.14|l . 

and H3.7|) . These transformation matrices are also unique corresponding to each 
Brunovsky form. 

3.4. Discussions. In the preceding subsection, we finished computing the two 
types of Brunovsky forms and the corresponding transformation matrices of a linearly 
controllable quadratic system, H1.5|l . Besides, we showed that the Brunovsky form 
of each type and the corresponding transformation matrices are unique. However, 
the uniqueness is dependent on r. That is, for different values of r, a system can 
have multiple solutions of type I or type II Brunovsky forms. A straightforward 
strategy to prevent the multiplicity is to set r = 0, and we have a simplified version of 
transformations, defined by (|1.8() . I.e., with (|1.8|l . a quadratic system has unique type 
I and type II Brunovsky forms. Further, one can follow the arguments of Kang and 
Krener [3, but with r = ^ and prove that H1.8(l still defines quadratic equivalence in 
the same sense. 

4. Computation of discrete quadratic Brunovsky forms. In this section, 
we apply the method proposed in the preceding section and solve Brunovsky forms for 
a discrete quadratic control system (|1.9() , but with simplified transformations defined 
in ll.8|l . We carry out our study in the same three-step, constructive procedure as 
for the continuous system. Compared to continuous system, the discrete system has 
one more term, and the relationships between coefficients and transformations and 
the resulted Brunovsky forms arc fundamentally different from those of continuous 
systems. 

4.1. Quadratically equivalent systems of ()1.9|l . Definition 4.1. Assume 
A g is in the linear Brunovsky form defined in {1-4^ , we define a linear operator 

L : SR"^" 5R"^" by 

L"P = P, 

LP = A^PA, 

p+ip ^ L-VP, i = 0,l,---. 

Properties of L. The linear operator L has the following properties: 

1. VP = when i > n. 

2. The nullity of L is 2n-l; if P G fcer(L), then Pij = when {n-i){n-j) > 
0, and elements in nth row and nth column can be arbitrarily selected. 

3. L is not invertible. 

4. If P is symmetric, LP is symmetric; but the inverse proposition may not 
be true. 

Remarks. Note that the linear operator L for the discrete system is different from 
that for continuous systems. This fundamental difference leads to the dif- 
ference in the relationships between coefficients and transformations and the 
Brunovsky forms. 

Theorem 4.2. Discrete system ^.fJ\} is quadratically equivalent, in the sense of 
under the quadratic transformations in to a system, whose ith (i — 1, - ■ ■ ,n) 



^Due to the difference in notations, r = is equivalent to saying /jl-'^l = in 
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equation is 



(4.1) x^{t + 1) = x^+i{t) + b^v{t) + x^{t)F,x{t) + Gix{t)v{t)+0{x, vf, 

where Xn+i{t) =0 is a dummy state variable, Fi, - ■ ■ ,Fn are symmetric, Gi is the ith 
row of the matrix G, and the coefficients F^ and hi (i — 1, ■ ■ ■ ,n) and G are defined 
by 



(4.2) 
(4.3) 
(4.4) 



Fi — Fi + Pi+i — LPi~biQ, 
G,^G,~ 2b'^P,A, 

hi — {Pi^nn, 



where Pn+i is a zero dummy matrix, bi is the ith element ofb, and (Pi)„„ is the {n, n) 
entry of the matrix Pi . 

Proof. The proof is similar to that of Theorem 13.21 and . therefore, omitted. ■ 
Remarks. Note that H4.3|l is different from H3.4|l . and we have an additional equation 
(El. 

4.2. A necessary condition for the equivalent systems. Lemma 4.3. The 
mapping from the coefficients G and G to the transformation matrices Pi (i = 1, ■ ■ ■ ,n) 
is given by 



(4.5) 



[P2 



[Pn 



(n) 



A= -G--G. 
2 2 



Proof. The proof is similar to that of Lemma Inland, therefore, omitted. ■ 
Lemma 4.4. The mapping from the coefficients Fi, Fi (i = 1, ■ ■ ■ ,n) and h to the 
transformation matrices Pi (i — 1, ■ ■ ■ ,n) and Q is given by 



(4.6) 

(4.7) 

and 
(4.8) 



3=0 j=0 

n-1 

Q = 5^£^"(F„_,-F„_,), 

j=0 



{Pl)nn = hi. 



, n - 1, 



i = !,■■■, n-1, 



Proof. The derivation of H4.6|l and H4.7|l from (|4.2() is similar to that in the proof of 
Lemma f3. 41 and, therefore, omitted. Here we only show how l|4.8|l can be derived as 
follows. From H4.6|l . we have (i = 1, ■ • • , n — 1) 



iPi+l)nn — (Ij*^1 )nn — '^^i^'' Fi-j)r. 

j=0 



(Pi) 



n — in — i ^ ^ (L'^^F^- 
3=0 
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Comparing this with 14.4|l . we then obtain (|4.8|l . from which we can see that the 
diagonal entries of Pi are uniquely determined by the coefficients Fi (i = 1, ■ • ■ ,n) 
and h. ■ 
Definition 4.5. With A e 3?"^" and L given by il.4\ ) and Definition \4-l\ ' 



spectively, we define a series of linear operators Xi ; 3?" 
by 



fi = o,i,---; 



(4.9) 



XoP = 



X,P 



{A^yxoP 



Properties of X^. The linear operators Xi (z = 0, 1, • • •) have the following proper- 
ties: 

1. Xi = when i > n. 

2. The rank of Xq is n(n + l)/2; denote the {i,j)th entry of P by Pij 
(i, j = 1, • ■ • , n), then the (i, j)th entry of XqP can be written as 



(4.10) 



(XoP),, = 



Pn-i+i,j-i+i, j>i; 
0, j < i. 



3. From Property 2, Xq is not invertible. Note that its continuous coun- 
terpart in Definition 13. 51 is invertible. 
Theorem 4.6. The mapping from the coefficients Fi, G, Fi, and G fi — 1, ■ ■ ■ ,n) 
to Pi is given by 



(4.11) 



n — 1 -| n — 1 I 

XoPiA = X^F,A +7^G-Y, XF,A - -G 



Proof. The proof is similar to that of Theorem 13 . 61 and omitted. 
Remarks. According to ()4.10fl . we can find the (i, j)th entry of XgPA, 



(4.12) 



(XoPiA), 



(Pl)ri-i-|-l,i-i, j > i] 

j<i. 



I.e., XqPA, the left hand side of H4.11|l . is an upper-triangular matrix, which 
contains all the non-diagonal entries of Pi. Therefore, Pi (i = 1, • • • , n) and G 
of a quadratically equivalent system of (|1.9|l must satisfy that the right hand 
side of H4.11|l be also upper triangular. Thus, H4.11|l is a necessary condition 
for all equivalent systems of (|1.9|l . including the Brunovsky form. 

4.3. Computation of the Brunovsky form and the transformation ma- 
trices. Theorem 4.7. In the Brunovsky form of J^1.9\) . the coefficients are 



(4.13) 



F_i = 0, Vi = l, 
G = 2{L + D), 



where L and D are the strictly lower triangular part and the diagonal part o/^"^/ XiFiA 



\G. 
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Proof. The first two terms on the right hand side of H4.11|l can be written as 

n— 1 

(4.14) ^X,F,A+-G = i + i^ + C/, 

i=l 

where L,D,U are strictly lower triangular, diagonal, and strictly upper triangular 
matrices respectively. To ensure the right hand side of (|4.11|) to be an upper triangular 
matrix, we set J2i=i '^i^i^ + to be 

n—l ^ 

(4.15) ^X,^,A+-G = i + i5, 

1=1 

in which there are n{n + l)/2 non-zero elements at most, and we hence have 

(4.16) XqPiA = U. 

Due to the properties of X^, no simple i^^'s with n(n+ 1) /2 non-zero elements can 
satisfy (|4.15|l . We thus simply set the coefficients as in (I4.13() . and obtain the only 
type of Brunovsky form of 1)1.9(1 : 

i 

(4.17) Xi{t + 1) ^ Xi+i{t) +biiy{t) + ^GijXj{t)i^ + 0{x,iyf, i=l,---,n. 

■ 

Remarks. The Brunovsky form of discrete systems corresponds to type II form of 
continuous systems and has n{n -\- l)/2 bilinear terms. After finding the 
Brunovsky form of a discrete quadratic system, we can solve the correspond- 
ing transformation matrices Pi {i — 1, ■ ■ ■ ,n) and Q as follows. From H4.12|l 
and H4.16|l . we can find the non-diagonal elements of Pi, and the diagonal 
elements from ((4.8(1 . Then Pi {i — 2, ■ ■ ■ ,n) and Q can be calculated from 
lira and llOl . 

5. Algorithms and examples. In this section, we summarize the algorithms 
for computing Brunovsky forms and the corresponding transformations of both con- 
tinuous and discrete linearly controllable quadratic systems with a single input. To 
prevent multiplicity in Brunovsky forms, we use the simplified version of transforma- 
tions, ((1.8(1 . for both systems. Thus, r = for all formulas in Section 3. Following 
each algorithm, an example is given. 

Algorithm 5.1. Computation of continuous quadratic Brunovsky forms. 

Step 1 Compute X^^ {^Y^^=i ^iPi + '^"'^ decompose it into L + D + U . 

Step 2 Compute Pi = L + D + L'^ and Ai = Xq{U - L'^). 

Step 3 For type I Brunovsky forms, set G — Q and solve the equation 

i=l 

to compute Fi (i = 1, ■ ■ ■ ,n — 1) as in 1^3.15]) and L'i.16]} . Go to Step 5. 
Step 4 For type II Brunovsky forms, set Fi = (i — \, ■ ■ ■ ,n) and G = 2Ai. Go to 
step 5. 
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Step 5 Compute (i — 2, ■ ■ ■ ,n) and Q from b\) and j^. ?| ). 

Example 1 . Find type I Brunovsky form of a continuous quadratic control sys- 
tem, which is already in type II Brunovsky form, 



(5.1) 



6 

6 



6 + o(e,/if, 



Solutions. In this system, n = 2, and coefficients are 

Fi=F2^ 0, 




G = 



1 



(5.2) 



Following Algorithm we find type I Brunovsky form, 



Xi 



X2 



2-^2 



X2 = i^ + 0{x,iy)^, 
and the corresponding transformations. 



(5.3) 



Xi, 

v. 



2-^27 



Reversely, if given 1)5. 2|l . we find the transformations for obtaining (|5.1|l . 

xi = 6, 

a;2 = 6 - 2'?2: 
1/ = ^, 

which are the inverses of (|5.3|l . 

Algorithm 5.2. Computation of discrete quadratic Brunovsky form. 
Step 1 Compute X)"=i^ XiFiA + anc? decompose it into L + D + U . 
Step 2 Compute G^2{L + D). 

Step 3 Compute the non-diagonal elements of Pi from \4.1'^ and \4.16}) . 

Step 4 Compute the diagonal elements of Pi from i4.8{ ). 

Step 5 Compute Pi (i — 2, ■ ■ ■ ,n) and Q from from j-^.d] ) and j^.Tj ). 

Example 2. Find the Brunovsky form of the following discrete system: 

ut + 1) = 6(0 + eAt) + e2{t) + m'(o + o(e, 

6(i + l) = /i(t)+M'(i) + 0(e,Ai)'- 
Solutions. In this system, n=2 and coefficients are 



Fi 

F2 = 0, 
G = 0, 



1 
1 
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Following Algorithm 15.21 we find the Brunovsky form, 



xi{t + l) = X2{t)+0{x,iy)^, 
X2it + 1) = iyit) + 0{x,L')\ 



and the corresponding transformations 



aw 



xi{t) + 2xl{t)+xl{t), 

X2it)-xlit)+xl{t), 

iy{t)-xl{t). 



Thus, this system is linearized. 



6. Conclusions. In this paper, we proposed a method for computing Brunovsky 
forms of both continuous and discrete quadratic control systems, which are linearly 
controllable with a single input. Our approach is constructive in nature and com- 
putes Brunovsky forms explicitly in a three-step manner, with the introduction of 
linear operators L and (i = 0, • • • , n — 1): (i) we first derived relationships between 
the coefficients of original systems, those of the quadratically equivalent systems, and 
the corresponding transformation matrices; (ii) after simplifying these relationships, 
we obtained a necessary condition for all equivalent systems; (iii) from the necessary 
conditions, we finally computed the Brunovsky forms and the corresponding trans- 
formations. However, the linear operators L and and, therefore, the Brunovsky 
forms are fundamentally different for continuous and discrete systems. For a con- 
tinuous quadratic system, which has at most n{n + 1)^/2 nonlinear terms, there are 
two types of Brunovsky forms with n{n — l)/2 nonlinear terms each; for a discrete 
quadratic system, there is only one type, and the numbers of nonlinear terms of a 
general system and its Brunovsky form are n{n+l)'^ /2+n and n(n-|-l)/2, respectively. 

For continuous systems, we used the same transformations as in |3] and found 
that the inclusion of transformation vector r yields multiple solutions of a type of 
Brunovsky forms. Therefore, we suggested to set r = in order to maintain the 
uniqueness. I.e., under transformations defined by (|1.8|l for both continuous and dis- 
crete systems, the Brunovsky forms and the corresponding quadratic transformations 
always exist and are uniquely determined by the original systems. In this sense, 
our study can be viewed as a constructive proof of the existence and uniqueness of 
Brunovsky forms. 

Finally, the method proposed in this paper could be extended for studying quadratic 
control systems, not linearly controllable or with multiple input. In addition, the 
method could be useful for applications involving analysis and control of quadratic 
systems. 
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